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*0 ■ Abstract 

O 

Q\ ■ We develop a formalism for evaluation of the transverse momentum dependence of cross 

sections of the radiation processes in medium. The analysis is based on the light-cone 

path integral approach to the induced radiation. The results are applicable in both QED 

and QCD. 



It is well known that at high energies the multiple scattering can considerably modify 
cross sections of the radiation processes in medium [l], [|. Recently this effect (called 
the Landau-Pomeranchuk-Migdal (LPM) effect) in QED and QCD has attracted much 
attention ||, f|, ||, |, [7|, || ||, |l(| (see also JTTJ] and references therein). In we have 
developed a new rigorous light-cone path integral approach to the LPM effect. There 
we have discussed the p^-integrated spectra. For many problems it is highly desirable to 
have also a formalism for the p^-dependence of the radiation rate. In the present paper 
we derive the corresponding formulas. Similarly to || our results are applicable in both 
QED and QCD. 

For simplicity we describe the formalism for an induced a — > be transition in QED for 
scalar particles with an interaction Lagrangian L int = \[wi)ci>a + i)£i>ci>b] (it is assumed 
that m a < rrib + m c , and the decay a — > be in vacuum is absent). The S'-matrix element 
for the a — > be transition in an external potential reads 

(bc\S\a) =ijdtdvXi)l{t,v)i)* c {t,v)i) a {t,T), (1) 

where ipi are the wavefunctions (incoming for i = a and outgoing for % = b,c). We write 

ipi as 

V'ifor) = -7if|fexp[-i(t-*)pi,J&(t,r). (2) 

\/ihji 

We consider the case when the particle a approaches the target from infinity, and normalize 
the flux to unity (it corresponds to \<fi,j\ — 1) at z — — oo for i = a and at z = oo for 
i — b,c. The case when the particle a is produced in a hard reaction in a medium (or 
at finite distance from a medium) will be discussed later. At high energy, E t ^> m i5 
the dependence of (pi on the variable r = (t + z)/2 at t — z =const is governed by the 
two-dimensional Schrodinger equation 

dcj)i 
i-Q- = H i<fii , (3) 

*=-U& + « u+ £.' (4) 

where /ij = p i>2 , p is the transverse coordinate, e^ is the electric charge, and U is the 
potential of the target. 

In the high energy limit from (|TJ), (0) one can obtain for the inclusive cross section 



d n 



a 



dxdc{ b dc{ c (2tt) 



-Rejdp 1 dp 2 j dz 1 dz 2 g(F(z 1 ,p 1 )F*(z 2 ,p 2 )), (5) 



Z\ <Z2 



where F(z, p) = 4>l(t, r)(f)*{t, r)<fi a (t, r)\ t=z , q bc are the transverse momenta, x = Pb, z /p a ,z 
(note that for the particle c p Cj2 = (1 — x)p a ,z), 9 = A 2 /[167rx(l — x)E%], (...) means 
averaging over the states of the target. Since the wavefunctions enter ([5]) only at t = z, 4>i 
can be regarded as functions of z, and p. In the Schrodinger equation @ z will play the 
role of time. We represent the z-dependence of 0« in terms of the Green's function, Ki, 
of the Hamiltonian (H). Then, diagrammatically, (0) is described by the graph of Fig. la. 



We depict K t (K*) by — > (<— ). The dotted line shows the transverse density matrices at 
large longitudinal distances in front of (z = Zj) and behind (z = Zf) the target .[] We will 




(b) 



Figure 1: The diagram representation of the inclusive spectrum (5) (a), and (7) (b). 

consider first the q c -integrated spectrum. For the sake of generality we assume that all the 
particles are charged in this case. Later we will give the formula for the totally inclusive 
spectrum when at least one of the final particles has a zero charge. For the q c -integrated 
case the transverse density matrix for the final particle c is given by a ^-function, and 
taking advantage of the relation 



Jdp 2 K(p 2 ,z 2 \p 1) z 1 )K*(p 2l z 2 \p 1) z l ) =5{p x - p[) 



(6) 



one can transform the graph of Fig. la into that of Fig. lb. The corresponding analytical 
expression reads 



d 3 



a 



z f z f 



;Rejdp b jdp' b jdp b dp' b dp a dp' a dp aii dp' a>l exp[-iq 6 (p 6J - p bJ )]jdz 1 jdz 2 



dxdq b (2n) 2 

X gSb(p b J, PbJ, Z f\Pb, Pb, Z2)M(p b , p' b , Z 2 \p a , p' a , Zi)S a {p a , p' a , Z 2 \pa tl , p a>i , Zi) , (7) 

where 

Si(p 2 , p 2 , z 2 \p 1: p[, z x ) = (K i (p 2 ,Z2\pi,z 1 )K*(p 2 ,z 2 \p' 1 ,z 1 )) (8) 

is the evolution operator for the transverse density matrix, and the factor M is given by 

M(p 2 ,p' 2 ,z 2 \p 1 ,p' 1 ,zi) = (K b (p 2 ,z 2 \pi,z 1 )K c (p' 2 ,z 2 \p 1 ,z 1 )K*(p 2 ,z 2 \p 1 ,zi)). (9) 

We assume that the target density does not depend on p. Then a considerable part 
of calculations can be done analytically. In (||), (§J) we represent the Green's functions 
in the path integral form. In the corresponding path integral formulas for Si and M 
the interaction of the particles with the target potential after averaging over the target 
states turns out to be transformed into the interaction between trajectories described 
by the Glauber absorption factors. For Si the corresponding absorption cross section 
is given by the dipole cross section an of interaction with the medium constituent of ii 
pair. The absorption factor for M involves the three-body cross section a a b c depending 
on the relative transverse vectors T bc = p b — p c and r ab = p a — p b . The factor Si can 
be evaluated analytically. The corresponding formulas are given in [ 12 , || . The factor M 



1 Strictly speaking, in (p, (J5|) the adiabatically vanishing at \z\ ~ \z% f\ coupling should be used. For 
simplicity we do not indicate the coordinate dependence of the coupling. 



after the analytical path integration over the center-of-mass coordinates can be expressed 
through the Green's function K bc describing the relative motion of the particles b and c 
in a fictitious abc system. The formula for M can be obtained from that given in || by 
replacing the final transverse coordinate p 2 by p 2 f° r the particle c. The expression for 
the probability of the a — > be transition at a given impact parameter which we obtain 
integrating analytically over all the transverse coordinates (except for r& = p b j — p b j) in 
(0) has the form 



d 3 P 



ReldTb exp(-iq b T b )jdz 1 jdz 2 g$ b {T b , z 2 )K bc {T bl z 2 \0, zi)$ n (T a , z x ) , (10) 



dxdq b {2ix) 2 ^J 
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where 




§a{T a ,Zi) =exp 


(Taa{T a ) I J ( „\ 


/ a«n\«) 




Zi 



<5>b(T b ,Z 2 ) 



exp 



°~bb( T b 



z f 



dzn(z) 



;n) 



are the eikonal initial- and final-state absorption factors,^ T a = xr b . 
for the Green's function K bc reads 



H, 



1 







in\Z aabc(T bc , T ab 



be 



2/j bc \dr 



be. 



The Hamiltonian 



(12) 



where n bc = E a x(l - x), r ab = -[r a + (1 - x)T bc ], L f = 2E a x(l - x)/[m b (l - x) + 
m 2 c x — m 2 a x(l — x)] is the so-called formation length. In (|TT|), (0) n(z) is the number 



density of the target. If the target occupies the region < z < L one can drop in (|10|) 
the contribution from the configurations with z\^ 2 < and Z\p > L. This follows from 

L 6c 



the relation for the Green's function K?„ for the Hamiltonian (1127) in vacuum 



l^bcT 

2z 



z 

T 



f\ 



o 



oo oo 

Re fdzK^ c (r, z\0, 0) = Re fdz (iy^-) exp 



This allows one to rewrite (0) in another form 

d 3 P 2 r 

dx^ b =j2^ Re J dTbeXp{ - iqbTb) 
z f 
dzi dz 2 g{$b(Tb,z 2 ) K bc (T b ,z 2 \0, z^ - K bc (r b , z 2 \0, zi) <b a {T a ,zx) 



(13) 



x 



z\ 



+ [M-rb, z 2 ) - 1] K° c (T b , z 2 \0, Zl ) [$ a (r a , zi) - 1]} 



(14) 



which demonstrates explicitly that the configurations with z\ y2 < and z Xy2 > L do not 
contribute to the radiation rate. Equations flT0|), ( |I4"D establish the theoretical basis for 
evaluation of the transverse momentum dependence of the LPM effect. 

2 Note that appearance of the eikonal absorption factors in (fLOl) is a nontrivial consequence of the 
specific form of the evolution operators S a< b |13], and is not connected with applicability of the eikonal 
approximation in itself. 



The integration over q 6 in ( |T4] ) gives the x-spectrum 



dP 

dx 



z f z f 



2Re dzi dz 2 g K 6c (0, z 2 \0, z x ) - # 6 ° c (0, z 2 \0,z 1 ) 



(15) 



which we obtained earlier in pi. There it has been derived using the unitarity connection 
between the probability of the a — > be transition and the radiative correction to the a — > a 
transition. The latter is described by the diagram of Fig. 2a which in turn using ([]) can 
be transformed into the graph of Fig. 2b. It can also be obtained directly from the graph 
of Fig. lb after integrating over q 6 . 





(b) 



Figure 2: The diagram representation of the radiative correction to the probability of 
a — > a transition. 



In the low density limit (14) can be written through the light-cone wavefunction \E^ C 

d 3 P 

dxdq b (2n 



l — Jdrdr' exp(-zq 6 r')^ c *(x, r - r')Ta- abc (r bc , r ab )^ b a %x, r) , (16) 



where T bc = r, T ab = —[(1 — x)t + xt], T = J dzn(z). This formula can be obtained 
from ([14]) taking advantage of the representation for \l/^ c through K bc 



*JT(*,t) = tA [ dzKl(r,0\0,z) 

AE„\ hx(l - x) -L 



established in [Q. Being divided by T (jig) gives a convenient formula for the Bethe-Heitler 
cross section in terms of the light-cone wavefunction. It worth noting that (|T6"D (and (|T0D, 



(0) as well) is valid if one can neglect the transverse motion effects on the scale of the 
medium constituent size. This assumes that the typical value of \z 2 — Z\\ in (|T0|), (|T4"D, 
which can be regarded as the formation length associated with the a — > be transition at 



3 The diagram of Fig. 2b gives only the term ex Kbc in ( |15|) (the corresponding integral is divergent in 
itself). Nonetheless, it yields the same result as (|l5|). Indeed, by adding and subtracting the contribution 
from the configurations z\ < Zf < Z2 , one can represent the contribution of the vacuum term as a sum of 
the imaginary term connected with the radiative correction to m a (which oc {zf — z{)) and the real term 
related to the wavefunction rcnormalization. The latter comes from the configurations Z\ < Zt < z 2 - 
This boundary effect is absent if the coupling vanishes at large \z\. In this case the vacuum term in 
( |l5| ) does not affect the x-spectrum. It is, however, convenient to keep the vacuum term to simplify the 
troublesome z-integration in (fL5|). Again, it allows one to use a constant coupling. 



a given q 6 , L/(g&), is much larger than the size of the medium constituent. If the LPM 
effect is not very strong the Lf(q b ) can be estimated replacing m bc by m 2 c + q^ in the 
above formula for Lf. Note that for Lf(q b ) 3> L the radiation rate can be written through 
\l/* c for arbitrary target density. In this case the target acts as a single scattering center 
and flTlD can be written in a form like fllTf ) but with the product Ta n bc being replaced by 



2{1 



exp 



\To, 



aba 



L 



This representation generalizes the formula for the pr-integrated 
spectrum derived in [|13 |. 

In general case one can estimate the radiation rate using the parametrizations a^ = 

— (Tab + t&c))- Then the Hamiltonian 



^ii T i > &abc = ^ab T ab + ^bc T bc + C ca T ca (here T ca = — 

(|12|) takes the oscillator form with the frequency Q(z) 

C ab (l — x) 2 + C bc + C ca x 2 . The Green's function for the osci 

z-dependent frequency can be written in the form 



(i-O 

V2 



n(z)C{x) 



1/2 



where C(x) 



E a x(l-x) 

lator Hamiltonian with the 



K n 



Ta,^Ti,«i 



7(^1^2) 
2ni 



expji a(zi,z 2 )rl + /3(z 1j z 2 )t 2 1 -^(z 1 ,z 2 )t 1 t 2 } . (17) 



The functions a, f3 and 7 in (JlTf) can be evaluated in the approach of Ref. [0. Then we 
can integrate analytically over T& in (ITO) , and represent the radiation rate as 



d 3 P 1 f i(z u z 2 ) 

axaq b \^) J Q{zi,Z2j 



exp 



i<4 



i{z\ - z 2 ) 



4Q(z u z 2 ) 



where the factor Q(zi, z 2 ) can be expressed through the parameters Cy, functions a, /3, 
7, and n. The formula for this factor is cumbersome to be presented here. 

Consider now the case when the particle a is produced in a medium or at finite distance 
from a target. Equation ( |10D holds in this case as well but now z% equals the coordinate 
of the production point. Given the representation ([10|) taking advantage of (|13D one can 
obtain a formula similar to (0) but with [$ a — 1] being replaced by $ a in the second term. 
Note, however, that, due to infinite time required for the formation of \l/^ c , equation fll6|) 
(and its analogue for arbitrary density at Lf(q b ) ^> L) does not hold in this case. 

Let us discuss briefly the totally inclusive radiation rate. It can be evaluated almost in 
the same way as the q c -integrated spectrum if one of the final particles has a zero charge, 
as this occurs for the e — > je transition in QED. Consider the case when e c = 0. Since the 
particle c does not interact with the medium the graph of Fig. la can be transformed into 
a graph like that of Fig. lb but with the propagator K c being connected to the lower abc 
vertex through the density matrix of the particle c. The corresponding formula (which is 
the analogue of flTU|)) reads 



d 5 P 



dxdq b dq c (27r) 2 



Re / dT h dr c exp [-i(q b r b + q c r c )] 



z f z f 



x jdzi dz 2 g®b(r b , z 2 )K hc {r h - t c , z 2 \0, zi)<& (t , zi) 



(19) 



where r a = xr b + (1 — x)t c . The ^-integration in ( [Of) can also be written as in fll4]) . In 
the low density limit and at Lf(q b ) ^> L the initial- and final-state interaction vanish. For 



this reason the analogue of ( |I~6"D and a similar equation for arbitrary density at £/(<#,) 3> L 
which can be obtained from fll~9|) are valid even when all the particles are charged. 

The generalization of the above results to the realistic QED and QCD Lagrangians 
reduces to trivial replacements of the two- and three-body cross sections, and vertex 
factor g. The latter, due to spin effects in the vertex a — > be, becomes an operator. The 
corresponding formulas are given in || |i~5|1 . 

The formalism developed can be applied to many problems. In particular, in QCD 
this approach can be used for evaluation of high-p^ hadron spectra, the p^-dependence 
of Drell-Yan pairs and heavy quarks production in ZiA-collisions, angular dependence of 
the parton energy loss in hot QCD matter produced in AA-collisions. It is also of interest 
for study the initial condition for quark-gluon plasma in AA-collisions. Some of these 
problems will be discussed in further publications. 
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